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Abstract. We study the problem of ensemble equivalence in spin systems with 
short-range interactions under the existence of a first-order phase transition. The 
spherical model with nonlinear nearest-neighbour interactions is solved exactly both 
' for canonical and microcanonical ensembles. The result reveals apparent ensemble 

inequivalence at the first-order transition point in the sense that the microcanonical 
entropy is non-concave as a function of the energy and consequently the specific heat is 
negative. In order to resolve the paradox, we show that an unconventional saddle point 
^ , should be chosen in the microcanonical calculation that represents a phase separation. 

O I The XY model with non- linear interactions is also studied by microcanonical Monte 

Carlo simulations in two dimensions to see how this model behaves in comparison with 
, the spherical model. 
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1. Introduction 

In statistical mechanics, we prepare an ensemble of macroscopic systems and calculate 
thermodynamic quantities by taking the average over the ensemble. When the system 
is isolated from the environment, the total energy is kept constant and the principle 
of equal weights leads to the microcanonical ensemble. On the other hand, when we 
consider a heat bath attached to the system to allow an energy exchange, we have 
the canonical ensemble characterized by temperature. These ensembles are generally 
considered equivalent and their thermodynamic potentials are related by the Legendre 
transformation. 

Equivalence of ensembles has been proven rigorously for systems with short-range 
interactions [1]. For systems with long-range interactions, there is no guarantee that 
two ensembles produce the same results in the thermodynamic limit. Typical examples 
include gravitational systems [2]-[Z] and fully-connected mean- field spin models [8]- 
[20] . For a review, see [21]. In the latter models, in particular, the interplay of long- 
range interactions and first-order phase transitions is now known to lead to ensemble 
inequivalence, typically as negative specific heat in the microcanonical ensemble. 

In systems with short-range interactions, by contrast, ensembles are equivalent in 
the thermodynamic limit and there should exist no anomalous effects except in finite- 
size systems [22]. In the present paper, we solve the multi-component spin model with 
nonlinear interactions in two and three dimensions exactly for the spherical model and 
numerically for the XY model. These models have been known to have first-order phase 
transitions in two and three dimensions [23]- [27]. We show that ensemble equivalence 
should be taken with special caution in these systems. 

The organization of this paper is as follows. In section 121 we define the model. The 
spherical (many-component) limit is solved exactly in section |3l The results for the 
canonical and microcanonical ensembles are compared. To study the system with finite 
component spins we use microcanonical Monte Carlo simulations in section HI The last 
section is devoted to summary and conclusion. 

2. n- vector model with nonlinear interactions 

We study the generalized n- vector model (0(?T,)-symmetric model) 



on a d-dimensional hypercubic lattice. The spin variable Si at site i is an n-component 
vector with the constraint Sf = I]"=i(>S'j")^ = n. The sum in the Hamiltonian is taken 
over nearest-neighbour pairs. The number of spins is = L*^, where L is the linear 
size of the system. In the standard n- vector model with linear interactions, the function 
V{x) is equal to x. Here, following [211 [25], we consider the form 




(1) 




(2) 
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The linear interaction is recovered if we choose p = 1. 

The hnear model in the limit n — )■ cxd is the ordinary spherical model and can be 
solved exactly. We shall call the nonlinear model also the spherical model for simplicity. 
The canonical analysis of the linear case is found in standard textbooks [281 129] . The 
microcanonical analysis was performed in [30] and [31]. We generalize their calculations 
to the nonlinear case. 



3. Spherical limit 

In the canonical ensemble, the generalized n-vector model ([1]) can be solved exactly 
in the spherical limit n oo [2^. The problem is reduced to solving a saddle-point 
equation for auxiliary variables. We solve the nonlinear model in the microcanonical 
ensemble and compare the results with those for the canonical ensemble. 



3.1. Saddle point equations 



(3) 



First, we briefly review how the problem is solved in the canonical ensemble 
following [25]. The partition function is written as 

r r In 

Z = Tr^exp pjnY,V{Si- Sj/n) Y[6[S 

where /3 is the inverse temperature and the trace denotes integrations over the spin 
variables. In order to carry out the integrations, the 6 function is expressed by a Fourier 
integral over the auxiliary variable Zi. We also introduce two kinds of variables pij 
(= Si ■ Sj/n) and \ij (to impose the constraint pij = Si ■ Sj/n) and write 

N 

Y[ dzi Y[ dXijdpij Tr exp 



i=l 



/3JnJ2V{p.,)-j:^^iSl 

{ij) i 



n] 



J2 ^ij ^^Pij 
{ij) 



Si ■ S i 



N 



W dzi W dXijdpij exp 

{ij) 



i=l 



PJn^V {pij) + n^Zi-nY^ XijPij 

{ij) i {ij) 



(4) 



+n In Tr exp —'^ZiS^ + '^ XijSiS 

\ i {ij) 

The spin trace is just a Gaussian integral over unconstrained scalars {Si} and can be 
evaluated using the lattice Green function [28l EH]. In the limit n — )■ oo, auxiliary 
variables are determined from the saddle-point equations. Following the conventional 
procedure used for the spherical model with linear interactions, we neglect the subscript 
dependence of the variables, X = Xi, z = Zi, p = pij (yi,j). Then, we can write 



exp 



n 

Nnd(3JV (p) + Nnz - NndXp + :t In G{k, z) 

^ k 
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where z = z/X. The lattice Green function in the momentum space is given by 

Gik,~z) = —^ rr- (6) 

z - T,t,=i cos kf. 

If we take the thermodynamic hmit, the sum over k is replaced by an integral as 

From the expression ([5]), we determine the state of the system by a set of saddle- 
point equations, 

X = /3JV'ip), 2\ = gCz), dp = z-^, (8) 



9C^) = I 7^ . .... • (9) 



where 



The auxiliary variables A and p are eliminated to obtain 

(3J=^^ ^. (10) 

z 



2V' - 



d dg{z) J 

For a given /3J, z is determined from this equation. Then, the free energy density 
/ = F/Nn is given by 



z 

d dg{l 

1 



+ -(l + ln27r). (11) 

Thus, by solving the simple saddle-point equation ( ITOl) for z, we can calculate the free 
energy as a function of {3. 

Next, we derive the equations in the microcanonical ensemble. If we compare (fTlll 
with the relation — /3F = —I3E + S", we may guess that the energy and entropy are given 

as 

s = -i 1,1 + i / ^ In G(*;, i) + i(l + In 271), (13) 

where e = E/NnJ and s = S/Nn. Here, 5 is obtained as a function of e by f|T2|) 
to determine the entropy s = s{e). These expressions ( |T2|) and ( |T3|) for the energy 
and entropy can indeed be derived directly from the microcanonical number of states 
Q = e^ = Tt 6{E ~ H) using the integral representation of the delta function. Following 
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the same procedure as in the canonical case, we find 

dt 



2tt 



exp 



iNnJt {e + dV{p)} + Nnz - NndXp 



N 



\ i=i {ij) J. 

Then, we impose the saddle-point conditions for the auxiliary variables to obtain the 



above result (|T2|) and f|T3|) as described in more detail in [Appendix A 



We are ready to study the ensemble dependence of system properties by comparing 
the canonical result (ITO!) . ( ITTj) and the microcanonical (fT2|) . (1131) . In the following, we 
focus ourselves on the cases of = 2 and d = 3. 



3.2. d = 2 




Figure 1. Saddle-point equation 
([TO)) in the canonical ensemble at 
d = 2. The function diverges 
logarithmically at the origin. 



Figure 2. Saddle-point equation 
p2p in the microcanonical ensem- 
ble at d = 2. The function is finite 
at the origin. 



Let us first write the lattice Green function ([6]) as, using the variable = 2(z — d), 

G(k, z) = 2 . (15) 

m2-f E^(2sin(V2))^ 

This is a decreasing function of m and the value at the origin m = determines 
the infrared behaviour. When d = 2, g{z) defined in iQ diverges logarithmically at 

m = {z = 2). 

We plot the right-hand sides of (ITO!) and (fT2|) in figures [1] and [21 respectively, 
for several values of p. In the canonical case, the right-hand side of (JTOi) diverges at 
z = d = 2 [m = 0) and is monotonically decreasing when p < 4. Therefore, in this 
case, for a given /3, z is determined uniquely. When p > 4, in a certain range of /3, 
three solutions exist and m is not determined uniquely. This is understood as the 
emergence of a first-order transition [25]. On the other hand, the function in (fT2|) for 
the microcanonical ensemble is finite at the origin (m = 0) and is a decreasing function 
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for arbitrary p. Since the functional value at the origin corresponds to the ground-state 
energy, the solution is determined uniquely for a given e larger than or equal to the 
ground-state energy. Nothing singular happens in this case. 




/(T) in the canonical ensemble for ^^le thermodynamically irrelevant 

d — 2, p ~ 1. (unstable or metastable) saddle- 

point solutions. 




Figure 5. Entropy density s(e) 
and inverse temperature /3(e) = 
ds{e)/de in the microcanonical 
ensemble for d — 2, p — 1. 




Figure 6. s(e) and /3(e) in the 
microcanonical ensemble for d = 
2, p — 5. s(e) is non-concave 
and correspondingly /3(e) is non- 
monotonic. The dotted parts 
correspond to the metastable and 
unstable branches in the canonical 
ensemble in figure ID 



From the obtained saddle-point solution in the canonical ensemble, we plot the 
free-energy density / for p = 1 and 5 in figures [3] and HI respectively. We see that 
a first-order transition occurs when p = 5. We also plot the entropy density s and 
the inverse temperature /3 = ds/de in the microcanonical ensemble for p = 1 and 5 in 
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figures E] and El respectively. For p = 5 in figure El we see that (3 is non-monotonic 
(the entropy is non-concave) for —0.21 < e < —0.07 and consequently the specific 
heat is negative. In this sense, ensembles may seem inequivalent. This behaviour is 
similar to the case of systems with long-range interactions, where the mean-field picture 
applies. A remarkable fact is that this mean-field-like behaviour has been found by 
exact calculations for the two-dimensional system with short-range interactions. 



3.3. Phase separation 

Ensemble equivalence is recovered in the present system if we choose a proper saddle- 
point solution which represents a phase-separated state. We have assumed in f|T3|) that 
the auxiliary variables are independent of the subscripts i and j. It implies that the 
phase is uniform in space. In order to describe the situation with phase separation, we 
divide the system into two parts with A^^i and N2 spins, respectively. The particular 
shape of the two sub- regions is not important, as far as they are geometrically compact 
objects, with a surface-to-volume ratio that vanishes in the thermodynamic limit (for 
instance, a cubic lattice may be divided into two slabs). We set auxiliary variables in 
each subsystem as z^^\ p^^\ X^^'> and z^'^\ p^'^\ X^'^\ respectively. In the thermodynamic 
limit, we expect that the interface terms between two subsystems are irrelevant due to 



the short-range nature of the system. We prove it rigorously in Appendix B Then, the 



number of states Q is written as the sum of contributions from two subsystems as 




+ nlnTrexpj - ^ z^'\sf^f + ^ \^'^ sf^ Sf 



(16) 



where Sf'^ {Sf''') represents the spin variable in subsystem 1 (2). The saddle-point 
equations are written as 

e = ^e, + -e„ (17) 

where zi = z^^^/X^^^ and Z2 = z^'^'^/\^'^\ Then, the entropy density is expressed as 

s=—s^ + —S2, (19) 

s. - -^ln.(.«) - ^/ l^ln - E cosO + \ + 1^2.. (20) 
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Figure 7. Comparison of tlie entropies obtained from the uniform and phase- 
separated solutions. The entropy from the phase-separated solution is represented 
by the straight lines, (a) When the entropy for the uniform solution is concave, the 
phase-separated solution is irrelevant because it has a smaller value of the entropy 
sjj than the true Sa- (b) When the entropy is non-concave, a first-order transition 
occurs and a phase-separated state is realized between ei and 62 because it has a 
larger entropy. 



A notable fact is that these expressions hold for any ratio of the separated phases, 
c = Ni/N, 1 — c = N2/N as well as for any choice of ei,2 and si,2- Thus, we should 
discuss what values of these parameters are actually chosen for a given fixed value of the 
energy e = ea- First, if the total entropy is concave, the hypothetical phase separation 
means that the value of the entropy would be s'^ = csi -|- (1 — c)s2, which is lower than 
the true entropy Sa as can be understood from figure [7](a). Thus, there is no phase 
separation in the stable state. Technically, this means that the exponent of the integral 
for the number of states Q = becomes largest at the saddle point representing the 
uniform state, not at the point corresponding to the phase-separated state. 

On the other hand, the situation is different when the total entropy is non-concave 
as we show in figure [Tl^b). At e = ea, the state with the entropy s = Sa = csi + (1 — c)s2 
is more stable than that with s'^ and is reahzed as the phase-separated state in the usual 
sense. Technically, the saddle point corresponding to this former state has the largest 
contribution to the integral. Thus, we can obtain the phase-separated state by relaxing 
the uniformity condition of the saddle-point solution. It should be noticed that only the 
microcanonical solution needs this non-uniform prescription of the saddle-point values. 
The uniform solution for the canonical case ([8]) shows no inconsistencies. 



Let us next consider the three-dimensional system, in which case g{z) is finite at 
m = {g{z = 3) ~ 0.505) and is monotonically decreasing. As shown in figures M 
and [U the saddle-point equation has no solution at low temperature or low energy. To 
avoid this difficulty, the zero mode = in ([9]) should be separated from the integral. 



3.4. d = 3 
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Figure 8. Saddle-point equation 
(fTU)) in tlie canonical ensemble at 
d = 3. 




10 go 



Figure 10. Zero-mode part of 
the saddle-point equation in the 
canonical ensemble at d = 3. 




Figure 9. Saddle-point equation 
(|12p in the microcanonical ensem- 
ble at d = 3. The value at the ori- 
gin is not equal to the ground state 
energy. 




Figure 11. Zero- mode part of 
the saddle-point equation in the 
microcanonical ensemble at d = 3. 



similarly to the Bose-Einstein condensation, as 

9(5)^,(* = 0.5)+9(5) = ^— i + Zp^JTlJ:;^. (21) 

The parameter z approaches 3 in the thermodynamic limit so that the first term gives 
a finite contribution go in this limit. Then, we can find the solution of the saddle-point 
equation by the replacement 

g(z)^go + gi3). (22) 

As depicted in figures [10] and [11], g^ can be fixed from the saddle-point equation for a 
given /3 (or e) below (or above) the values achievable in figure [H] (or figure [9|). Hence, 
there exists a solution for any /3 or e, the latter being larger than or equal to the 
ground-state energy. 



We plot the free energy in the canonical ensemble in figures [12] and [13] for p = 1 
and p = 5, respectively. A continuous transition to the zero-mode condensed phase 
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Figure 12. /(T) in the canonical 
ensemble for d = 3, p = 1. The dot 
denotes the transition point, where 
go starts to be finite. 




Figure 13. /(T) in the canonical 
ensemble for d — 3, p — 5. 




Figure 14. s(e) and /3(e) in the 
microcanonical ensemble for d = 
3, p = 1. The dot denotes the 
transition point. 




Figure 15. s(e) and /3(e) in the 
microcanonical ensemble for d = 
3, p — 5. The putative continuous 
transition point denoted by the 
dot is hidden by a first-order 
transition. The dotted parts 
correspond to those in figure 1131 



is observed for p = 1. It is replaced by a discontinuous transition for p = 5. From 
the microcanonical analysis, we plot s and /3 in figures [H] and [15] for p = 1 and 
p = 5, respectively. As can be understood from these figures, transitions between 
the condensed and non-condensed phases exist in both ensembles in three dimensions. 
It is discontinuous for p = 5. Similarly to the d = 2 case, the apparent ensemble 
inequivalence (negative specific heat in the microcanonical ensemble) can be avoided by 
the proper choice of the saddle-point solution. 
4. Monte Carlo analysis 

In order to check if the results of the previous section is specific to the spherical model, we 
study the n = 2 model, the XY model, in two dimensions by Monte Carlo simulations. 
A canonical Monte Carlo analysis of the present model has already been carried out 
in [23], and the microcanonical case has been done in [33]. However, these previous 
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studies are not sufficient to clarify the problem of ensemble equivalence, and we analyze 
the same model here from our own point of view. We refer the reader to |3l] for a 
related study, where the ensemble equivalence was examined in a mean-field model on 
random graphs. 

The original XY model with linear interactions in two dimensions exhibits the 
Kosterlitz-Thouless transition [35] . In |36], a detailed study is reported on how the 
transition can be changed in nonlinearly-interacting systems. Since our aim is not to 
study this topological transition, only the caloric curve has been calculated in our Monte 
Carlo analysis. 




Figure 16. /3(e) from the micro- 
canonical Monte Carlo calculation 
for n = 2, c? = 2 and p = 1. The er- 
ror bars are smaller than the sym- 
bol size. 



40 

3 

30 
20 
10 




L=20 
L=40 



-0.10 -0.08 -0.06 -0.04 8 -0.02 0.00 

Figure 17. /3(e) for n = 2, d = 2 
and p = 40. The error bars are 
smaller than the symbol size. 
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Figure 18. /3(e) for n = 2, d 2, 
p = 40, and L = 50 near the non- 
monotonic region. We have taken 
the average over 20 independent 
runs. 




-0.08 -0.07 



-0.06 ^ -0.05 



-0.04 -0.03 



Figure 19. Energy histogram 
in the phase coexistence region 
from the canonical Monte Carlo 
calculation. We have chosen n — 2, 
d ^ 2, p = 40, L ^ 50, and 
/3 = 19.35. 
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Several system sizes have been analyzed, L = 20, 40, 50 and larger in some cases. 
To study the microcanonical ensemble, we exploit the demon algorithm by Creutz |32j . 
In this algorithm, for a given energy e, the demon energy Ed is calculated so that the 
sum of the system and demon energies is kept constant. Then, the inverse temperature is 
obtained from the probability distribution Prob(i?D) ~ exp(— /3ii^D)- We have performed 
10^ Monte Carlo steps per spin for each run. 

In figures [16] and [T71 we plot the energy dependence of the inverse temperature 
for p = 1 and p = 40, respectively. For p = 1, we see that /3 is a monotonically 
decreasing function of e. It is different for p = 40, where the function shows a non- 
monotonic behaviour. Figure HH] highlights this property for L = 50. For a given /3, e 
is not determined uniquely in a narrow region, which suggests the existence of negative 
specific heat. We have performed canonical Monte Carlo calculations using the simple 
Metropolis algorithm to see the energy histogram, and the result is depicted in figure [191 
which clearly shows that a first-order transition exists in the form of phase coexistence. 
We have confirmed that the non-monotonic region of the caloric curve remains up to 
the size L = 100 in the microcanonical simulations. An extrapolation suggests that 
it would persist to the thermodynamic limit. Thus, the negative specific heat seems 
to exist in the microcanonical ensemble also in the two-component system as in the 
spherical model. 

We speculate that a reason for the apparent ensemble inequivalence for p = 40 
in Monte Carlo simulations may be that a phase separation, as discussed above for 
the spherical model, has not been realized in our simulations because of a very long 
relaxation time: The system has to spontaneously break up into two spatially separated 
regions with different macroscopic states, which could take a very long time to be realized 
in the microcanonical simulations. 

5. Summary and conclusion 

We have studied the n- vector model (O(n)-symmetric model) with nonlinear short- 
range interactions in two and three dimensions. The exact solution of the spherical 
model shows ostensible inequivalence of canonical and microcanonical ensembles through 
negative specific heat in the latter ensemble. We have argued that this paradox can 
be resolved by explicitly taking into account a phase separation, which increases the 
entropy (thus increases the thermodynamic stability) in the microcanonical ensemble. 
It is noticed that the proper choice of the saddle-point solution is required in the 
microcanonical ensemble to represent the state with phase separation. Such a solution 
must be considered when the uniform ansatz of the saddle-point solution gives a non- 
concave entropy. We note that this procedure is needed only in the microcanonical 
ensemble. In the canonical ensemble, the uniform solution is sufficient to represent the 
stable state of the system. Another interesting aspect is that the exact solution of the 
spherical model is similar to the mean-field solutions applicable to long-range interacting 
models in the sense that the calculations of the (canonical) partition function and 
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the (microcanonical) entropy reduced to steepest descent integrations. The difference 
between the short- and long-range systems is apparent when we consider the phase- 
separated state. The interface term cannot be neglected in the long-range system, which 
can be understood as an important source of ensemble inequivalence. Our analysis has 
succeeded to highlight this difference in an exactly solvable example. 

The XY model with nonlinear interactions has been shown to behave similarly by 
Monte Carlo simulations in two dimensions. Discrepancies between ensembles may in 
this case be due to the long relaxation time to the fully-stable phase-separated state in 
the microcanonical simulation. We expect that a more elaborate method such as the 
one developed in [37] may resolve this problem. 

The ra-vector models with nonlinear short-range interactions have been known to 
have the unusual property of the existence of first-order phase transitions even in two 
(and higher) dimensions [23]- [27]. We have identified an additional highly non-trivial 
property of apparent ensemble inequivalence, which we expect to stimulate further 
studies of this very unusual class of models. 

Appendix A. Derivation of the entropy in the microcanonical ensemble 

We consider the number of states for a given energy E 

Q = Ti6{E-H) = J ^Tre*(^-^)*. (A.l) 

This expression has a similar form to the partition function in the canonical ensemble 
except for the integral over t and the factor e*"^*. We may thus replace P in the partition 
function by it. Therefore, the calculation goes along the same line as in the canonical 
case and we can write 

^= f ^ f Y[dziY[dXijdpij exp it| E Jn^ V (pij) j + n ^ 
i Hi) \ (ii) J i=l 



/ TV 

KjPij + n In Tr exp - ^ ZiS^ + ^ XijSiSj 



(A.2) 

Then, we impose the uniform ansatz for Zi, pij and Xij and obtain the number of states 



-n 

(ij) \ «=1 (ij) 



as 



Q = exp 



it(^E + NndJV{p)) + Nnz - NndXp 



(A.3) 



/ N 

-|- 77. In Tr exp — ^ zSf + ^ XSiSj 
\ i=i (ij) 

and the saddle-point conditions 

e + dV{p) = 0, X = ttJV'{p), 2X = g{z), dp = z - (A.4) 
where e = E/NnJ. Combining these results, we finally obtain f|T2|) and ([1 



Ensemble equivalence in spin systems with short-range interactions 



14 



Appendix B. Bounds on interface effects in the free energy of the Gaussian 
model 

In order to rigorously justify the calculations using only the two phase-separated regions 
without interface terms, let us estimate the order of magnitude of the effects that the 
interface terms have on the free energy of the Gaussian model. We define the free energy 
of the Gaussian model as 

/ (1) (2) \ 

h{x) = In Tr exp [-zJ^Sf + XJ2 ^i^J + ^^Y. ^i^:! ' (^-1) 

where the summation with superscript (1) runs over all interactions within the two 
independent (phase-separated) subsystems and the summation with superscript (2) is 
for interactions across the interface. The parameter A will be assumed to be positive 
without losing generality on a bipartite lattice. Notice that we have assumed that the 
interactions have common values in the two subsystems. We will show later that this 
restriction can be removed. The boundary conditions are assumed to be free in the x 
direction and periodic otherwise. Here the term 'interface' stands for the region in the 
middle of the system that runs perpendicular to the x axis and separates two subsystems, 
whereas the 'boundary' is for the outmost sites of the total system. Equation ( IB.ip 
indicates that the interface interactions have the strength Xx and all other interactions 
have A. 

Our goal is to prove that 

\h{l) - h{0)\ < cN,, (B.2) 

where A^f, is the number of interactions across the interface and c is a quantity 
asymptotically independent of Nf, and N (total number of sites). This inequality (IB. 20 
shows that the presence and absence of boundary interactions affect the free energy only 
by a term proportional to and thus can be neglected in the thermodynamic limit 
where the leading term is of order A^. 

Let us first notice that the derivative of h{x) is non-negative, the first Griffiths 
inequality, 

h'i-) = A > 0, (B.3) 

where —H{x) is the effective Hamiltonian appearing in the exponent of (IB.ip . The 
denominator of f IB.SP is positive. The numerator is also non-negative for < x < 1: 



Each term of the expansion of the numerator 

(2) oo /(I) (2) \" 

^rY.S^S,Y.-T{T.S^S,+xY.S.S^ e-^^^'f (B.4) 

{ij) «=0 (ij) J 

is composed of integrals of the form 

TrS^Sp'^e-'^'^l (B.5) 
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which is zero (if any one of a,b,c,-- - is odd) or positive (otherwise). The second 
derivative is also non-negative: 

2. / . /o^ . \ 21 



h"{x) = A' 




' / (2) ^ 



> 0, 



where (■ ■ ■)g stands for the average by the weight e Thus h'{x] 

and is bounded by h'{l){> 0) for < x < 1. Therefore 



\h{i)-hm 



/•I , dh 






dh 


/ dx— 


< 


/ dx 


Jo dx 




Jo 


dx 



< h'il). 



(B.6) 

is non-decreasing 
(B.7) 



Our task is to upper-bound h'{l). From the definition of h{x), this derivative is expressed 
as 

h\l) = XNk{SiS,)G{x = 1), (B.8) 

where (ij) is a bond across the interface. According to (lB.7p . if we are able to prove 
that {SiSj)G{x = 1) is finite in the thermodynamic limit (A^ — oo,Nb — oo), we will 
have finished the proof that h{l) and h{0) are no more different than a quantity of order 
A*";,. This implies that the contribution of the interface interactions can be neglected in 
the computation of the bulk free energy. 

Finiteness of r{y = 0) = {SiSj)G{x = 1) can be shown as follows. The argument y of 
r{y) stands for the strength of interactions connecting the left-most sites and right-most 
sites along the x direction. In other words, y = 1 corresponds to the periodic boundary 
and y = is for free boundary in the x direction (Remember that x = 1 ensures that 
the interactions across the interface exist). All other directions have periodic boundary 
conditions. The Hamiltonian is modified as 

(1) (2) (3) 

-H{x = l,y) = -zY: Sfj + S^S, + A ^ S,S, + yA J] S,S„ (B.9) 

* (ij) (ij) (ij) 

where the final sum with superscript (3) runs over the boundary bonds. Let us assume 
for the moment that we have proved the following inequality, 

<r{y = 0) <r{y = 1). (B.IO) 

Since r(l) is the single-bond correlation for fully-periodic boundary conditions, we can 
calculate it explicitly by taking the derivative of the free energy with respect to A and 
diving the result by the total number of bonds. The explicit form is available for this 
quantity in ([5]) and it is easy to see that r(l) is positive and finite provided that z > d. 
This ends the proof that r(0) is finite. 

To prove ( IB. 101) . we first notice r(0) > 0, the first Griffiths inequality, which can 
be proved as we did above. Next we take the derivative of r{y), 



(3) 

r'{y) = ^Yl {{SiSjSiSm)G 

(Im) 



{SiSj)G{SlSm) 



G) 



(B.ll) 
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The definition of (■ ■ ■)g is sliglitly modified in that the Hamiltonian —H{x = l,y) is 
now used. Since the integral defining (■ ■ ■)g is Gaussian, Wick's theorem apphes, 

{SiSjSiSm)G — {SiSj)G{SlSm)G 

= {S,Si)g{SjS^)g + {S,S^)g{SjSi)g > 0, (B.12) 

due to the first Griffiths inequahty. The proof of ( IB. 101) thus completes. 

Finally, we show that the result apphes also to the case where the two subsystems 
have different values of A. Let us replace A by mA (0 < m < 1) for one of the two 
subsystems. The other subsystem keeps the original value of A. Then, r is a function 
of y and u. The derivative of r{y, u) with respect to u has an expression very similar 
to (IB. Ill) , which can be shown to be positive as before. Thus, r{y,u) < r{y,l) for 
< u < 1. Since u = 1 is for the system already treated above, we know that r(l, 1) 
is finite. It then follows that r{0,u){< r(0, 1) < r(l, 1)) is finite. All other parts of the 
proof can trivially be generalized to accommodate < m < 1. Q.E.D. 

The condition of the outer boundary (free or periodic) along the x axis can also be 
shown to be irrelevant in the thermodynamic limit. To outline the process, let us define 

= lnTre-^(^=^'^). (B.13) 

The goal is to prove 

|J(1)-J(0)| <ciV,, (B.14) 

where c is a quantity that converges to a finite value in the thermodynamic limit and 
A'";, is the number of bonds appearing in the summation with superscript (3). To show 
this, according to our experience above, we should prove the relations 

f{l) = cm, j\y)>0. (B.15) 

These can be proved in the same manner as before. 
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